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Abstract 

The magnetic entropy change [ASm(T, HU around the phase transition temperature Tọ is in- 
vestigated by the scaling method for Feo.5Coo.5Si, which exhibits a skyrmion phase below Te. The 
parameters of ASy(T, H) exhibit field dependent behaviors. The ASm(T, H) curves under high 
field can be well scaled into a single universal curve independent of external field and temperature. 
However, ASm(T, H) curves under low field become divergent just below Tc, which indicates a 
characteristic of first-order transition. The scaling investigation of ASm(T, H) curves indicates 
that the phase transition in Feo.5Coo.5Si is of a weak first-order type in low field region, while it is 
driven into a second-order one under high field. This weak first-order phase transition in low field 
region resembles that in typical skyrmion system MnSi which is caused by the critical fluctuation. 
The result suggests that critical fluctuation plays an important role in the phase transition and 


formation of skyrmion state. 
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I. INTRODUCTION 


The chiral magnetic materials have attracted much attention recently due to the exotic 
physical phenomena, such as magnetic quantum phase transition, anomalous Hall effect, 
non-Fermi liquid behavior, and tricritical point [1-4]. These fascinating phenomena origi- 
nate from the competition and balance between the spin, orbital, charge, and lattice degrees 
of freedom. In B20 compounds belonging to the P213 space group, due to the noncentrosym- 
metric structure, the lack of inversion symmetry usually results in the Dzyaloshinsky-Moriya 
(DM) interaction, which is 1~2 orders of magnitude smaller than the ferromagnetic coupling 
[5]. The competition between DM spin-orbital coupling and strong ferromagnetic exchange 
usually generates a helical magnetic ordering with long-wave length modulation [6-8]. The 
helical magnetic ground state can be polarized into conical ordering by the external field [9]. 
When an appropriate external magnetic field is applied, a skyrmion state, which is a topo- 
logically stable vortex-like spin texture with nano-size, appears at the boundaries between 
the conical ordering and paramagnetic phases just below To [9, 10]. The B20 compounds, 
such as Fe;_,Co,Si, FeGe, MnSi with the cubic cell, are typical systems to investigate the 
chirality of the spin helix and the quantum phase transitions [11-14]. 

In the chiral magnetic systems, Fe;_,Co,Si has been prominent for its easy modulation 
through the doping of Co [14, 15]. As we know, FeSi is paramagnetic with a strongly 
correlated narrow gap while CoSi is a diamagnetic metal [16]. However, Fe;_,Co,Si ex- 
hibits itinerant magnetic properties [17, 18]. In addition, the magnetic helical ordering in 
Fe,_,Co,Si can be changed from the right-handed to left-handed tuned by the content of 
Co [19]. The spin-interaction can also be modulated by the doping of Co [20]. It has 
been demonstrated that the change of magnetic chirality and spin interaction by doping are 
due to the tune of the DM interaction [15]. Especially, the half-doped Fey. sCoo.,;Si exhibits 
complex non-collinear magnetic ordering and precursor phenomena [13, 21]. Therefore, the 
investigation of Fep.5Cog.5Si is of great importance to uncover the magnetic interaction of 
this system. 

As we know, the magnetic entropy changes when the phase transition occurs, which is an 
effective means to study the magnetic interaction [22-25]. In this work, we perform a scaling 
investigation of the magnetic entropy change around the phase transition temperature for 


Fep,.;Cog.55i. The scaling investigation of the magnetic entropy change suggests that the 


phase transition in Fep.;Cog5Si exhibits a weak first-order characteristic under low field, 


which is driven into a second-order one by high magnetic field. 


II. EXPERIMENT 


A single crystal sample of Hoas Coa säi was prepared by the Czochralski method [26]. The 
physical properties were checked in elsewhere [20]. The magnetization measurement was 
performed by using a Quantum Design Vibrating Sample Magnetometer (SQUID-VSM). 
The initial isothermal magnetization around the phase transition temperature (Tc) were 
recorded to calculate the magnetic entropy change. Before the data collection, the magnetic 
field was swept to zero in an oscillating mode. Then, the temperature was warmed to 100 
K (T > Te) waiting for two minutes. Then the temperature is decreased to the target 
temperature. After waiting for two minutes, the no-overshoot mode was applied to get 
a precise magnetic field. For the measurement of each initial magnetization curve, the 
magnetic field is decreased to zero field in the oscillating mode to remove the remained field. 
Each initial isothermal magnetization curve was measured in the same procedure to ensure 


a precise magnetic field and temperature. 


III. RESULTS AND DISCUSSION 


Figure 1 (a) shows the temperature dependence of magnetization |M(T)] (left axis) for 
Feo .5Cog,5Si under zero-field-cooling (ZFC) and field-cooling (FC) sequences with an applied 
field H = 10 Oe. The magnetic transition temperature Je is determined from the peak 
of M(T) curve as To ~ 44 K. The temperature dependence of reciprocal susceptibility 
[x '(T) = (M/H)"!] curve is also given in Fig. 1 (a) (right axis). The y~'(T) curve 
exhibits a straight linear behavior above Tc, which indicates a paramagnetic behavior. The 
x~1(T) obeys the Curie-Weiss law: x(T) = C/(T — Bear). where C = 0.0467(9) [emu K/g 
Oe] is the Curie constant, and dcw = 49.12(6) K is the Curie-Weiss temperature. The inset 
of Fig. 1 (a) plots the isothermal magnetization [M(H)] at 4 K, which displays a typical 
magnetic ordering behavior with no magnetic hysteresis. In order to generate the magnetic 
entropy change, the initial M(H) curves in the vicinity of Tc are given in Figure 1 (b). The 
inset of Fig. 1 (b) gives the M(H) curve in the low field region at 20 K, which shows that 


saturate magnetic field Hs ~ 500 Oe. 
Based on the thermodynamical theory, the magnetic entropy change LA Gu, H)] induced 
by the external field is given by: [27]: 
ASu(T,H) = IK ES GEN (1) 
with the Maxwell’s relation [0S(T, H)/OH]r = [OM(T, H)/OT|n, ASu(T, H) is written as 
[28]: 


Aut Deere / AH 2) 


where H is the external magnetic field, T is the temperature, and H™ is the maximum of 
external magnetic field. Based on the initial M(H) curves in Fig. 1 (b), the temperature 
dependence of ASy [—AS),(T)] curves under different H are calculated. Figure 2 (a) 
and (b) show the —ASj,(T) curves around Te under high and low field, respectively. For 
each -ASy(T) curve under selected field, a peak appears around Te. The position of the 
peak shifts slightly to the higher temperature direction with the increase of H, with the 
temperature offset AT ~ 2 K. The value of -ASy(T) almost decreases symmetrically on 
both sides. The magnitude of each -ASy(T) curve for the fixed temperature increases with 
the applied magnetic field and reaches at the maximum value of 0.74 J/(kg-K) when H = 30 
kOe. 

The shape of —ASm(T) curve changes as the increase of field, such as the broadening 
and increase of the curves. The change of -ASy(T) vs. H can be quantified by the 
parameters, which usually follows power laws as a function of field. Figure 3 (a), (b), and 
(c) show the field dependence of parameters of |AS77°”| (the maximum of the magnetic 
entropy change), d6rwxam (the full width at half maximum), and RCP(S) (the relative 
cooling power), respectively. It can be seen that all parameters of |AS7i°*|, Aru: par, and 


RC'P(S) are field dependent, which obey the relations as [29]: 
Age oc H” 
Öörwum X H : (3) 
RCP(S) x H° 
where n, b, and c are exponents related to the magnetism of the system. For |ASY7“| vs. 


H, it can be obtained that n = 0.7431(2). Previous investigations on different soft magnetic 


amorphous alloys have demonstrated that n is close to 1 well below T¢ while it approaches 
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2 well above Teo [30]. However, n exhibits a complex behavior around Tc, which is related 
to the critical behavior. Based on the Landau mean field theory, n should be equal to 2/3 
at To. However, for a typical soft magnetic amorphous alloy, n ~ 0.75 around Te [30]. The 
obtained n = 0.7431(2) for Feo.5Coo.5Si is very close to that of soft magnetic amorphous 
alloy, which is also confirmed by no magnetic hysteresis phenomenon on M(H) curve. The 
Orwum implies the broadening of -ASy(T) curves with the increase of H, as shown in Fig. 
3 (b). The fitting of the rwym vs. H gives that b = 0.3817(6). The RC’ P(S) is defined as 
[30, 31]: 

RCP(S) = [ASK x Orwam| (4) 


The fitting of RCP(S) vs. H in Fig. 3 (c) yields that c = 1.1034(4). 
The variation of the phase transition can be also demonstrated by the magnetic specific 
heat change [AC p(T, H)]. The ACp(T, H) is obtained from the change of ASu(T, H) as 


[32]: 
OASM(T, H) 
ee: (5) 


The temperature dependence of —AC’p under different H is depicted in Fig. 4. The -ACp 


AC, (T, H) = OST, H) — C,(T,0) = 


changes from positive for T Te to negative for T < Te with temperature decreasing, 
which corresponds to paramagnetic to magnetic ordering phase transition. The magnetic 
transition temperature can be determined by temperature corresponding to the zero point of 
ACp [ACp(T, H)|r. = 0]. However, it can be seen that Te exhibits a fluctuation of AT ~ 2 
K, which is in agreement with the peak shift of ASm(T, H) in Fig. 2. 

The field dependence of magnetic entropy change |-ASy(H)] at different temperatures 
are depicted in Fig. 5 (a) and (b). The magnetic entropy change increases monotonically 
with the increase of external magnetic field. Figure 5 (a) and (b) depict the |-ASu(H)] 
above and below To, respectively. The maximum value of -ASy(T) appears around 46 K 
as shown in Fig. 5 (a). The value n can be calculated by [33]: 


din(ASu(T, Hu 


n(T, H) = din(H) (6) 


Figure 5 (c) and (d) show the temperature dependence of n under high and low field based 
on Fig. 5 (a) and (b). Then vs. T [n(T)] curve under high field is raised with the increase of 
field as shown in Fig. 5 (c), where n reaches the minimum value around To. However, n(T) 


curves under low field become divergent, as shown Fig. 5 (d). Generally, the divergence 


of the n(T) implies a first-order phase transition [34, 35]. The type of the phase transition 
should be firmly testified by the scaling of ASy(T, H) curves. 

The -ASy(T, H) curves can be investigated by the scaling method based on the 
critical phenomenon [29-31]. According to scaling laws, the experimental -ASu(T, H) 
should collapse onto two independent universal curves below and above Te with proper 
reference temperatures and normalization. The —ASm(T, H) is normalized by AS), = 
ASu(T)/ASY*(T)|, and the temperature is normalized by reference temperature T, into 
a re-scaled temperature 0. The horizontal coordinate is marked as the re-scaled temperature 


0, which is defined as [29]: 
T — Te 
~~ T, — To (7) 


By choosing proper T,, the equivalent points of the different experimental curves should col- 


H 


lapse onto the same point on the universal curve [29]. The experimental AS, should collapse 
onto two independent universal curves below and above To respectively. The combination of 
the two universal curves below and above Te just form a whole universal magnetic entropy 
change curve. After scaling and normalizing, T, should locate at a value of 6 = 1 on the 
universal curve [AS‘,(6)]. However, in multiple magnetic phases, the use of two reference 
temperatures is necessary. After normalizing the curves, T,ı and Ta are just corresponding 
to 0 = —1 and 0 = +1 respectively. The re-scaled temperature @ is defined as [29]: 
o 0- = (Te —T)/(Tn — Te), T < Te (8) 
6, =(T—Te)/(T2 — To), T > Te 
where 7}, and T, are the reference temperatures below and above Teo respectively, which 
corresponds to ASy(T,1, Tr2) = IAS. According to the scaling equation, the order of 
the magnetic phase transition can be judged by this scaling method. If the curves around Te 
can be re-scaled into a single universal curve independent of H and T, the phase transition 
is of a second-order type [34]. Otherwise, the divergence of the re-scaled A Gul) curves 
suggests a first-order one [22]. This method has already been extensively applied in many 
kinds of systems, such as MnSi, Gd(Si,Ge;_,)4, and many Heusler alloys [23, 36-38]. 
The construction of a universal scaling of the experimental -ASy(T, H) are plotted in 
Fig. 6. The ASu(T) curves for different H were normalized as A Su on the longitudinal 
coordinate, where the maximum of all curves are located on single point. The temperature 


is plotted as 0. Figure 6 (a) and (b) depict the AS7,(0) curves under high and low field, 
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respectively. In Fig. 6 (a), all experimental data collapse into a single universal curve, where 
T, and T, are field dependent as shown in the inset of Fig. 6 (a). The good convergence of 
the AS7,(0) curves under high field indicates the magnetic transition is of the second-order 
type [22]. A universal curve can be fitted, which is independent of the external measurement 
conditions but just determined by the intrinsic magnetization. The n values at high field 
can also be scaled into n(@) as shown in the inset of Fig. 5 (c), which confirms the scaling 
of -ASyu(T, H) curves at high field. However, in Fig. 6 (b), AS7,(0) curves below 900 Oe 
cannot be scaled into a single curve. The A Gulf) curves becomes divergent just below To. 
The non-convergent behavior of the scaling curves of magnetic entropy change suggests a 
first-order phase transition under low field below 9000e [22]. 

The scaling investigation of the magnetic entropy change in Fey 5Co9.,5Si demonstrates 
that the magnetic transition exhibits a weak first-order characteristic under low field, while 
displays a second-order one in high field region. In other word, the external field drives 
the weak first-order magnetic transition into a second one. In fact, the small angle neutron 
scattering (SANS) experiment has also suggests that the phase transition in Feo 5Coo.5Si 
exhibits a weak first-order characteristic under zero field [39], which is in agreement with 
the result in this work. This field-driven phase transition is usually observed in the B20 
compounds. In MnSi, which is a typical skyrmion system analogous to that in po son säi, 
the first-order phase transition is also be driven into a second one by the external field [22, 40]. 
According to the Brazovskii scenario, the first-order phase transition from paramagnetic to 
helimagnetic ordering in the B20 compound is induced by the critical fluctuation under zero 
field, which can be suppressed by the external magnetic field [41]. Therefore, it is suggested 
that the first-order characteristic induced by critical fluctuation plays an important role in 


the phase transition and formation of skyrmion in Fep,5Coo,.5Si. 


IV. CONCLUSION 


In summary, we have performed scaling investigation of the magnetic entropy change 
around Tc in the B20 compound Fe95C00,;Si. The parameters of ASm(T, H) curves 
(ASV, örwam, and RCP(S)) depends on the field as power law behaviors. The 
ASyu(T, H) curves under high field above 1 kOe can be well scaled into a single univer- 


sal curve AS7,(0) independent of external field and temperature. However, AS7,(0) curves 


under low field below 900 Oe become divergent just below Tc, which indicates a charac- 
teristic of first-order transition. The scaling investigation of the magnetic entropy change 
suggests that the phase transition in Feg.;Cog5Si is of a weak first-order type in low field 
region below 900 Oe, while it is driven into a second-order one under high field. This field- 
driven phase transition is analogous to that in MnSi, which is attributed to the critical 
fluctuation. The result suggests that the critical fluctuation plays an important role in the 


phase transition and formation of skyrmion in Fep,5Cop.5Si 
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FIG. 1: (Color online) (a) Temperature dependence of magnetization [M(T)] (left axis) and re- 
ciprocal susceptibility [x HTH (right axis) under H = 10 Oe for Feo,5Coo,5Si (the inset shows the 
magnetization as a function of field [M(H)] at 4K); (b) the initial isothermal M(H) around Tc 


(the inset gives the one at 20 K in the low field region). 


-AS, (J/kg.K) 
> > > 
N > Nn 


= 
= 


30 40 50 60 70 
T(K) 


FIG. 2: (Color online) (a) Temperature dependence of —ASjy |-ASm(T)] under different H in 
high field region; (b) —A Sm (T) curves under different H in low field region. 
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FIG. 3: (Color online) The field dependence of parameters of -ASy(T,H) curves: (a) -ASH 
vs. H, (b) örwum vs. H, and (c) RCP(S) vs. H (red solid curves are fitted). 
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FIG. 4: (Color online) Temperature dependence of —AC, [—AC,(T)] under different H for 
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FIG. 5: (Color online) (a) and (b): Field dependent of -ASu [-ASu(H)] at different tempera- 


tures above and below Te: (c) and (d): temperature dependence of n curves under high and low 


field (the inset of (c) shows the scaled n vs. 6 for the curves under high field). 
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FIG. 6: (Color online) Scaling of the -ASy(T, H) curves: (a) the normalized ASy(T,H) as a 
function of scaled @ in high field region (the inset gives Tu and Ta as a function of H); (b) AS7,,(0) 


curves in low field region (the inset depicts the magnified region below Te). 
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